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We present the results of extended theoretical LDA+DMFT calculations for a new iron-pnictide high 
temperature superconductor NaFeAs compared with the recent high quality angle-resolved photoemission 
(ARPES) experiments on this system [1]. The universal manifestation of correlation effects in iron-pnictides 
is narrowing of conducting bands near the Fermi level. Our calculations demonstrate that for NaFeAs the 
effective mass is renormalized on average by a factor of the order of 3, in good agreement with ARPES data. 
This is essentially due to correlation effects on Fe-3d orbitals only and no additional interactions with with any 
kind of Boson modes, as suggested in [1], are necessary to describe the experiment. Also we show that ARPES 
data taken at about 160 eV beam energy most probably corresponds to = tt Brillouin zone boundary, while 
ARPES data measured at about 80 eV beam energy rather represents fcz = 0. Contributions of different Fe-3d 
orbitals into spectral function map are also discussed. 


PACS: 71.20.-b, 71.27.+a, 71.28.+d, 74.70.-b 

1. INTRODUCTION 

The family of iron based high-temperature supercon¬ 
ductors first discovered in 2008 [2] still attracts a lot of 
scientific attention. Experimental and theoretical works 
on these materials are now discussed in several extended 
reviews [3, 4, 5, 6]. Detailed comparison of electronic 
band structures of iron pnictides and iron halcogenides, 
together with some related compounds was given in 
Refs. [7, 8]. 

One of the classes of iron pnictides is the so called 
111 system with parent compound Lii_a;FeAs with 
Tc=18 K [9, 10]. LDA band structure of the LiFeAs 
was first described in the Refs. [11, 12]. 

One of the most effective experimental techniques to 
probe electronic band structure of these and similar sys¬ 
tems is the angle-resolved photoemession spectroscopy 
(ARPES) [13]. A review of the present day status of 
ARPES results for iron based superconductors can be 
found in Ref. [14]. 

Soon after the discovery of iron based superconduc¬ 
tors it was shown both experimentally [15, 16, 17, 18, 19] 
(mainly by ARPES) and theoretically [20, 21, 22, 23] 
(within the LDA+DMFT hybrid computational scheme 
[24]) that electronic correlations on Fe sites are essen¬ 
tial to describe the physics in these materials. The main 
manifestation of correlations is simple narrowing (com¬ 
pression) of LDA bandwidth near the Fermi level by the 
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factor of the order from 2 to 4. At the same time the 
topology of ARPES determined Fermi surfaces is quite 
similar to those obtained from simple LDA calculations, 
showing two or three hole cylinders around P-point in 
the Brillouin zone and two electron Fermi surface sheets 
around (tt, tt) point. 

This work was inspired by recent high quality 
ARPES data for NaFeAs system [1] and is devoted 
to the detailed comparison of these results with 
LDA+DMFT calculations of electronic structure of this 
system, showing rather satisfactory agreement with 
these experiments. Thus, only the account of electronic 
correlations is sufficient to explain the major features 
of electronic spetrum of NaFeAs, and there is no need 
for any additional interactions with any kind of Boson 
modes (as was suggested in Ref. [1]). 

2. ELECTRONIC STRUCTURE 

The crystal structure of NaFeAs has tetragonal 
structure with the space group P4/nmm and lat¬ 
tice parameters a = 3.9494 A, c — 7.0396 A. The 
experimentally obtained crystallographic positions are 
the following Fe(2b) (0.75, 0.25, 0.0), Na(2c) (0.25, 
0.25, ZNa), As(2c) (0.25, 0.25, zas), zas=0.20278, 
ZiVa=0.64602 [25]. That is quite similar to LiFeAs crys¬ 
tal structure [9, 11]. 

In Fig. 1 we show LDA band dispersions (on the 
right) and densities of states (DOS) (on the left) cal¬ 
culated within FP-LAPW method [26]. Bands in the 
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Fig. 1. LDA calculated band dispersions (right) and 
densities of states (left) of paramagnetic NaFeAs. The 
Fermi level Ef is at zero energy. 


vicinity of the Fermi level have predominantly Fe-3d 
character and are essentially similar to the previously 
studied case of LiFeAs described elsewhere [11, 12]. The 
As-4p states belong to the -2 to -5 eV energy interval. 

To perform DMFT part of LDA+DMFT calcula¬ 
tions we used CT-QMC impurity solver [30, 31]. In 
order to link LDA and DMFT we exploited Fe-3d and 
As-4p projected Wannier functions LDA Hamiltonian 
for about 1500 /c-points. Standard wien2wannier inter¬ 
face [27] and wannierOO projecting technique [28] were 
applied to this end. The DMFT(CT-QMC) computa¬ 
tions were done at reciprocal temperature /3 = 40 with 
about lO’’ Monte-Carlo sweeps. Hubbard model interac¬ 
tion parameters were taken to be C/=3.5 eV and J=0.85 
eV as typical values for pnictides in general and close 
NaFeAs relative - LiFeAs in particular [32, 33, 34]. 

Fig. 2 shows the comparison of orbital resolved 
densities of states for Fe-3d shell of NaFeAs obtained 
within LDA (gray dashed line) and LDA+DMFT (solid 
gray and black lines). Solid gray and black lines show 
LDA+DMFT densities of states obtained by different 
methods of analytic continuation. Gray lines are ob¬ 
tained directly from DMFT(CT-QMC) Green function 
G(t) by maximum entropy method [29]. Overall line- 
shapes of LDA+DMFT densities of states are iden¬ 
tical to those already published in the literature for 
LiFeAs [32, 33, 34] and NaFeAs [34]. Most affected by 
correlations are Fe-3d(t2g) orbitals xy and degenerate 
xz,yz. These orbitals form narrow pronounced peaks 
near the Fermi level. On the other hand Fe-3d(eg) or¬ 
bitals 3^^ — and x^ — y^ just remind the broadened 
LDA densities of states. 

To produce LDA+DMFT spectral function maps for 
direct comparison with ARPES data we need to know 
the local self-energy E(a;). To find it we have to per¬ 
form analytic continuation from Matsubara frequencies 



Fig. 2. Comparison of orbital resolved densities of 
states for Fe-3d shell of NaFeAs obtained within LDA 
(gray dashed line) and LDA+DMFT (solid gray and 
black lines). Zero energy is the Fermi level. 

to real ones. To this end we have applied Fade approxi- 
mant algorithm [35]. The fact that both gray and solid 
lines coinside well in Fig. 2 tells us, that this analytic 
continuation is done rather satisfactory. Gorresponding 
self-energies for different Fe-3d orbitals near the Fermi 
level are shown on Fig. 3. From the real part of self¬ 
energy we can obtain the mass renormalization factor 
for different orbitals: m*/rrixy ~3.8, m*fmxz,yz ~3.9, 
TO*/m 3 j, 2_,.2 Ri2 and ra* jmx'^-y'^ ~1. These numbers 
agree well with variety of previous theoretical works for 
LiFeAs and NaFeAs [32, 33, 34]. Thus only the account 
of local Coulomb correlations on the Fe sites is enough 
to produce such renormalization and no extra interac¬ 
tion with possible Boson mode is necessary in contrast 
to the proposal of Ref. [1]. 

Typically experimental ARPES data are presented 
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Fig. 3. LDA+DMFT calculated self-energies for differ¬ 
ent Fe-3d orbils of NaFeAs near the Fermi level. Black 
lines - real part, gray lines - imaginary part. The Fermi 
level Ef is at zero energy. 
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Fig. 4. Comparison of experimental ARPES (left 
panel) [1] and LDA+DMFT (middle panel) spectral 
functions in the MFM high symmetry direction for 
NaFeAs for the wide range of binding energies contain¬ 
ing Fe-3d and As-4p states. On the right panel max¬ 
ima of experimental (crosses) [1] and theoretical (white 
lines) extracted from corresponding spectral functions 
are presented. The Fermi level Ef is at zero energy. 

in a rather narrow energy interval of few tenth of eV 
close to the Fermi level (for LiFeAs and NaFeAs see 
Refs. [34, 36, 37, 38, 39]). However, in Ref. [1] ARPES 
data were measured down to a quite large binding ener¬ 
gies about 6 eV with rather high resolution allowing to 
extract different bands. 

In Fig. 4 we compare experimental ARPES spectral 
functions for NaFeAs (left panel) [1] along the the MFM 
high symmetry direction with LDA+DMFT calculated 
(middle panel) spectral function map for a wide energy 
window. On both of these panels one can see rather 
high intensity region from 0 to 0.5 eV formed by quasi¬ 
particle bands near the Fermi level and then from -2 to 
-5 eV we can observe As-4p bands. To compare exper¬ 
imental and theoretical bands dispersions on the right 


panel of Fig. 4 we plot the dispersions for the maxima 
of experimental (crosses) and theoretical (white lines) 
spectral functions. 

According to Ref. [1] ARPES bands line shapes re¬ 
mind very much the LDA bands, compressed by an al¬ 
most constant factor of the order of 3 for all energies. By 
analyzing the real part of self-energies from Fig. 3 we 
can convince ourselves, that this correlation narrowing 
is essentially frequency dependent. Extended discus¬ 
sion of similar situation was given in our recent work on 
KaFe 2 Se 2 [42]. Actually, the LDA bands located in the 
interval from -0.5 eV to 0.25 eV become more narrowed 
due to correlations. At larger energies, the bands stay at 
about the same positions as in LDA or get more spread 
in energy since the slope of the real part of self-energy 
is changed to the positive one. 

As to As-4p bands ARPES experiment resolves only 
2 bands instead of 6 (2 As atoms in the unit cell). De¬ 
spite the general shape of the bands being quite similar 
in both cases, the experiment shows As-4p states about 
0.5 eV lower in energy than obtained in LDA+DMFT. 
This can be explained in the framework of generalized 
LDA’+DMFT calculations [40, 41], which allows one a 
better description of Fe(3d)-As(4p) energy splitting, as 
was shown for example for KaFe 2 Se 2 system [42]. In¬ 
deed our LDA’+DMFT calculations showed, that As-4p 
states appeared about 0.5 eV lower in energy. 

Between quasiparticle bands and As-4p bands there 
is a rather low intensity region (-0.5 eV - -2 eV) seen 
in Fig. 4 on the left and middle panels. First of all, 
it appears because there are almost no bands in this 
energy interval, and secondly in this region we have a 
crossover from the well defined quasiparticle bands with 
quite low damping to the rest of the bands placed at 
higher binding energies. This fact is illustrated by gray 
lines on Fig. 3, representing the imaginary parts E"(w) 
of LDA+DMFT calculated self-energies for all Fe-3d or¬ 
bitals. Near energy zero (Fermi level) E"(w) is about 
0.2 eV or less for all correlated states. At the same time 
real parts of the self-energies E'(u;) has negative slope 
near the Fermi level, which corresponds to well defined 
quasiparticles. Following S'(uj) behavior one can find 
that it has peak at about 0.25 eV, which corresponds to 
the end of quasiperticle region and E"(w) grows quite 
rapidly beyond this energy. Nearly the same behavior of 
i;'(w) and F"(uj) was assumed in the Ref. [1] and related 
to interaction with some “unknown Boson mode”, dis¬ 
tinguishing NaFeAs as unconventional superconductor. 
Again we claim that just the local Coulomb correlations 
on the Fe sites can do all that alone. 

Now we turn to the quasiparticle bands dispersions 
in the close proximity if the Fermi level. Correspond- 
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Fig. 5. Comparison of experimental ARPES (panels c and d) [1] and LDA+DMFT (panels b and e) spectral functions 
in the MFM and AZA high symmetry directions for NaFeAs near the Fermi level. On the panels a and f experimental 
(crosses) [1] and theoretical (white lines) maxima dispersions of spectral functions are presented. The Fermi level Ep is 
at zero energy. 


ing comparison of experimental ARPES data and the¬ 
oretical LDA-I-DMFT spectral functions for NaFeAs is 
given in Fig. 5. Here we show only the data for MPM 
high symmetry direction, the results for other symmetry 
drections can be found in the Supplementary Material 
[44]. Experimental data were obtained at two different 
rather distinctive beam energies around 80 eV and 160 
eV (see panels c and d on Fig. 5). At bird eye view for 
both beam energies experimental picture looks similar, 
but in fact there are some remarkable differences. For 
84 eV data xy and xz, yz bands close to the Fermi level 
are more intensive as compared to 160 eV data. On the 
other hand band at about -0.2 eV looks more 

intensive in 160 eV data. 

To clarify this fact we suggest following explanation. 
It is well known that by varying the beam energy in 
ARPES experiments one can access different values of 
kz component of the momenta [13]. However to get the 
precise value of kz one should know the exact geometry 
of ARPES experiment [13], work function and inner po¬ 
tential for this particular material [43]. Since we do not 
know all these precisely, we can try some speculations. 
In Fig. 5 we plotted LDA-l-DMFT calculated spectral 
functions for kz = d (panel b) and fcj, = tt (panel e). 
Now moving from fcz = 0 (panel b) to fcz = tt (panel 

e) we can observe the same trend as one goes from 84 
eV (panel c) to 160 eV (panel d) beam energy in the 
experiment. 

Although iron based superconductors have pro¬ 
nounced layered structure still these systems are quasi 
two-dimensional and thereby possess some finite disper¬ 
sion along kz axis. This fact is reflected on panels a and 
f of Fig. 5 where LDA-I-DMFT spectral function max¬ 
ima dispersions (white lines) are shown at fcz = 0 (panel 
a) and kz = tt (panel f). For the case of fcz = tt (panel 

f) xz, yz bands are no more degenerate in P-point. One 


of xz, yz bands branches goes down in energy to -0.2 eV 
and becomes degenerate with one of 3z^ — bands. At 
the same time 3z^ — band goes down to -0.4 eV at P- 
point and becomes more flat. All that, in contrast to the 
kz = 0 case, results in higher intensity of LDA-I-DMFT 
spectral function (panel e) around -0.2 eV and lower 
intensity at the Fermi level. The later one agrees bet¬ 
ter with 160 eV ARPES data, than with 84 eV ARPES 
data. Here one should stress that xy band and one of 
xz, yz bands branches right below the Fermi level keep 
their shapes almost unchanged for both = 0 and 

kz = TT. 

Somewhat larger intensity of LDA-I-DMFT spectral 
functions in comparison with experiment near M point 
arises because of quite strong xy contribution in this 
region. However in the ARPES data [1] xy band is al¬ 
most hidden, perhaps due to matrix elements effects (see 
Supplementary Material [44]). Note also that shown ex¬ 
perimental ARPES maxima (crosses on panels a and f) 
in accordance with Ref. [1] do not depend on beam 
energy. 


3. CONCLUSION 

In this paper we have presented the results of 
extended LDA-I-DMFT(CT-QMC) theoretical analy¬ 
sis of recent high quality angle-resolved photoemission 
(ARPES) experiments on a new iron-pnictide high tem¬ 
perature superconductor NaFeAs [1]. The well known 
and rather universal manifestation of correlation effects 
in iron-pnictides is the renormalization (narrowing) of 
conducting bands near the Fermi level by a factor of 2 
to 4. Corresponding mass renormalization factors for 
different orbitals were obtained from LDA-I-DMFT cal¬ 
culations and, in our opinion, no extra interaction with 
some “unknown Boson mode” distinguishing NaFeAs as 
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unconventional superconductor is necessary in contrast 
to the suggestion of Ref. [1]. 

Also we have shown that ARPES data taken at 160 
eV beam energy most probably corresponds to fcz = tt 
B rillouin zone boundary, while the data measured at 
about 80 eV beam energy reproduces = 0. Theoret¬ 
ical analysis of spectral weight redistribution support 
this point of view. Comparison of different Fe-3d or¬ 
bitals contributions to spectral function maps for verti¬ 
cally and horizontally polarized ARPES data also favors 
the last statement (see Supplementary Material [44]). 
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Supplemental Material to “Electronic structure 

of NaFeAs superconductor: LDA+DMFT 
calculations compared with ARPES 
experiment” 

In this Supplement we provide more results of 
our calculations for other symmetry directions in the 
Brillouin zone and some additional comparisons with 
ARPES experiments at different polarizations. 

In some sense MXM direction shown in Fig. 1 is 
the simplest one among others since only few bands are 
present here. The most intensive region here is around 
X point. For = tt 3z^ — band goes a bit down in 
energy (panel f). It leads to lowering of intensity around 
X point for theoretical spectral function (panel e) and 
quantitatively reproduce ARPES data at 160 eV (panel 

d). ^ 

In Fig. 2 for XPX high symmetry direction qualita¬ 
tive picture of bands evolution from kz = 0 to kz = tt is 
the same as for MPM direction (see Fig. 5 in the main 
text). Again the ARPES data at 80 eV (panel c) agrees 
better with kz = 0 LDA-bDMFT results (panels a,b). 
Most intensive spots of spectral function are formed at 
the crossing of xz,yz branches at -0.1 eV. While for 
kz = TT most intensive region appears around -0.2 eV, 
where and xz,yz bands are dominating. 

To discuss different Fe-3d orbitals contribution to 
spectral function maps we used experimental ARPES 
spectral functions obtained for different polarizations 
[1]. In Fig. 3 panel a corresponds to vertical polariza¬ 
tion ARPES data in the MPM high symmetry direction 
taken at 160 eV and panel f - to horizontal polariza¬ 
tion. For vertically polarized beam “cap”-like structure 
around P-point is formed mainly by xz, yz orbitals (pan¬ 
els b and d for LDA-I-DMFT results). Surprisingly the 
intensity of xy band (panels c and e) is quite low in 
ARPES data and even not addressed in Ref. [1]. The 
kz dispersion of these bands near the Fermi level is al¬ 
most absent. 

Horizontally polarized beam (panel f) wipes out 
3^2 _ ^2 forming “M”-like structure around -0.2 
eV. For kz = tt it has higher intensity than for kz = 0. 
It is in better agreement with to 159 eV data. 
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Fig. 1. Comparison of experimental ARPES (panels c and d) [1] and LDA+DMFT (panels b and e) spectral functions 
in the MXM high symmetry direction for NaFeAs near the Fermi level. On the panels a and f experimental (crosses) 
[1] and theoretical (white lines) maxima dispersions of spectral functions are presented. The Fermi level Ef is at zero 
energy. 



Fig. 2. Comparison of experimental ARPES (panels c and d) [1] and LDA+DMFT (panels b and e) spectral functions 
in the XFX high symmetry direction for NaFeAs near the Fermi level. On the panels a and e experimental (crosses) 
[1] and theoretical (white lines) maxima dispersions of spectral functions are presented. The Fermi level Ef is at zero 
energy. 
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Fig. 3. Comparison of experimental ARPES spectral functions with different polarization (panel a - vertical polariza¬ 
tion, panel f - horizontal polarization) [1] and LDA-I-DMFT spectral functions for different Fe-3d orbitals: panels b-e - 
xz,yz and xy contributions, panels g,h - 3z^ — contribution in the MFM high symmetry direction for NaFeAs near 
the Fermi level for = 0 and = tt cases. The Fermi level Ef is at zero energy. 


















